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Abstract 

The magnetic field probability P(B) is calculated from Ginzburg-Landau theory for various lattices of vortex lines 
in type-II superconductors: Ideal triangular lattices, lattices with various shear strains and with a super lattice of 
vacancies, and lattices of short vortices in films whose magnetic field "mushrooms" near the surface. 
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1. Introduction 

Type-II superconductors like Niobium and many 
alloys allow magnetic flux to penetrate in form of 
magnetic flux lines, i.e., vortices of the supercur- 
rent, each vortex carrying one quantum of magnetic 
flux. This effect was predicted in 1957 by Alexei 
Abrikosov [1], who got for this the Nobel Prize in 
Physics 2003. Abrikosov flux lines arrange to a more 
or less perfect triangular vortex lattice that exhibits 
interesting structural defects that may be calculated 
from Ginzburg-Landau (GL) theory [2] or by treat- 
ing the vortex lattice as a continuum with non-local 
elasticity [3] . The vortex lattice may even melt into 
a "vortex liquid" [4] . Pinning of vortices by material 
inhomogeneitics [5] together with thermal fluctua- 
tions of the vortices can cause a rich phase diagram 
in the magnetic field-temperature plane [6], with a 
melting line and an order-disorder line [7, 8] at which 
the weak elastic disorder ( "Bragg glass" ) suddenly 
changes to a plastically deformed or even amorphous 
vortex arrangement. The vortex lattice can be ob- 
served by decoration, magneto-optics, Hall probes, 
neutron scattering, magnetic force microscopy, and 
by muon spin rotation (/iSR). /iSR experiments can 



give valuable information about the vortex lattice, 
see the recent review [9]. 

When pinning and thermal fluctuations may be 
disregarded (e.g., in clean Niobium with very weak 
pinning) the vortex lattice exists when the applied 
magnetic field B a lies between the lower critical field 
B c i and the upper critical field B c i . In this interval 
the internal average induction B is smaller than B a , 
i.e., the magnetization M = B — B a is negative (dia- 
magnetic behavior), ranging from —B c \ to while 
B ranges from to B C 2- For B a < B c \ the super- 
conductor expels the magnetic field (B — 0, ideal 
Meissner state) and for B a > B C 2 the superconduc- 
tor is in the normal conducting state (B — B a ). This 
applies to long superconductor cylinders or slabs in 
parallel B a . For other geometries and for inhomo- 
geneous materials, demagnetization effects modify 
this picture and the magnetization curve in general 
has to be computed numerically, e.g., for thick or 
thin strips, disks, and platelets in perpendicular B a 
[10]. However, in the particular case when a homo- 
geneous specimen with the shape of an ellipsoid is 
put into a uniform applied field B a , then the de- 
magnetizing field (caused be the magnetization) in- 
side the ellipsoid is also uniform and superimposes 
to B a , thus generating an effective applied field 
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B, 



B a — 2VM(B.; N = 0). 



(1) 
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Solving Eq. (1) for the effective internal field Bj, one 
obtains M = M(B a ; N) = M(B 2 ; N = 0). Here N 
is the demagnetization factor, and the ideal magne- 
tization curve [M(B a ; 0) for N = 0, e.g., from GL 
theory] should be inserted. In general, TV is a ten- 
sor, but when B a is along one of the three principal 
axes of the ellipsoid, then N and all fields in (1) are 
scalars. For long cylinders or slabs in parallel B a one 
has N = 0, for spheres N = 1/3, for long cylinders 
in perpendicular field N = 1/2, and for thin plates 
and films one has 1 — N <C 1 . The sum of the N along 
the three axes is always N\ + N2 + N3 = 1 . For the 
Meissner state one finds B = 0, M(B a ;0) = -B a , 
Bi = B a /(l - N), and M{B a -N) = -B a /(l - N), 
which means the vortex penetration starts at the ef- 
fective penetration field B' cl = (1 — N)B C \ where 
Bi = B c \. In the interval B' cl < B a < B C 2 such a 
pin-free ellipsoid contains a perfect vortex lattice. 

In this paper the probability P(B) that at a ran- 
dom point inside the superconductor a muon sees an 
induction value B, is considered for ideal triangu- 
lar vortex lattices and for various possible perturba- 
tions of it, namely, various types of shear deforma- 
tion [11], a super lattice of vortex vacancies, the sur- 
faces of a film in a perpendicular field, and random 
displacements. This probability (or field density) is 
defined as the spatial average 



P(B') = (S(B'-B(r))) r . 



(2) 



In it B' is the independent variable, B(r) the spa- 
tially varying magnetic field, and 5(x) is the ID delta 
function, which for computation may be replaced by 
a narrow Gaussian whose width may depend on B' . 
One easily shows that 



J — i 



P(B)dB = 1, 



P(B)BdB = B, 



(3) 



2. 2D ideal triangular vortex lattice 

The local magnetic field B(x, y) of the vortex lat- 
tice may be calculated for all values of the reduced 
average induction b — B/B C 2 and Ginzburg-Landau 
parameter k by an elegant iteration method that 
minimizes the GL free energy F with respect to 
the Fourier coefficients of the periodic solutions 
B(x,y) and order parameter uj(x,y) = f 2 = \ip\ 2 
where tp(x,y) = f exp(iip) is the complex GL func- 
tion [12,13]. In the usual reduced units (length 
A, induction \/2B c , energy density B 2 /fi , where 
A is the magnetic penetration depth and B c = 
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Fig. 1. The field density P(B) of the ideal triangular vortex 
lattice obtained from GL theory for k = l/\/2 = 0.707 and 
b = 0.25 to 0.95. Here and in the figures below the variable 
B is normalized such that B = means -B m i n and B = 1 
means B m ax, and the inset enlarges the region near -B max - 




Fig. 2. The P(B) as in Fig. 1 but for for n = 1.4. 

B C 2/\/2k is the thermodynamic critical field with 
B C 2 — < J > o/(2tt^ 2 ), £ = X/k the coherence length 
and $o = h/2e = 2.07 • 10~ 15 Tm 2 the quantum 
of magnetic flux) the spatially averaged free energy 
density F of the GL theory referred to the Meissner 
state (tp = 1, B = 0) in the superconductor reads 



F _/ (i-M 2 ) a | 



v 



IK 



-A U 



B 



(4) 



Here B(r) = V x A, A(r) is the vector potential, 
and (...) = (1/V) J v d 3 r . . . means spatial averag- 
ing over the superconductor with volume V. Intro- 
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Fig. 3. The magnetic field variance a = ([B(x,y) — B] 2 ) 
of the triangular FLL for k = 0.85 to 200 plotted in units 
of B C 2 as yfa ■ (k 2 — 0.069)/S c 2 (solid lines) such that the 
curves for all k collapse near 6 = 1. The dashed lines show 
the same functions divided by (1 — 6) such that they tend 
to a finite constant 0.172 at 6 = 1. All curves are plotted 
versus Vb= ^B/B c2 to stretch them at small 6 values and 
show that they go to zero linearly. The upper frame 0.383 
is the usual London approximation. The limit for very small 
6 is shown as two dash-dotted straight lines for k = 5 and 
k = 10. The upper frame 0.383 shows the usual London 
approximation. 

during the super velocity Q(r) = A V<p/n and the 
magnitude /(r) = \ip\ one may write F func- 
tional of these real and gauge-invariant functions, 



(i-/ 2 ) 2 , (V/) : 



+ / 2 Q 2 + (VxQ) 2 ).(5) 



In the presence of vortices Q(r) has to be chosen 
such that V x Q has the appropriate singularities 
along the vortex cores, where / vanishes. By min- 
imizing this F with respect to ip, A or /, Q, one 
obtains the two GL equations together with the ap- 
propriate boundary conditions. For periodic lattices 
with one flux quantum per vortex, in the sense of a 
Ritz variational method one uses Fourier series for 
the periodic trial functions with a finite number of 
Fourier coefficients «k and 6k, 



;(r) = <Xk(1 — cosKr) 



K 

B(r)=B + b K cosKr, 

K^0 

Q(r) = Q4r)+ ^6 K ^^sinKr, 

K#0 



(6) 
(7) 

(8) 



where K = K mn = (K x ,K y ) are the reciprocal lat- 
tice vectors of the vortex lattice with positions R m „, 



R« 



(mil + nx 2 ; ny 2 ) , 



K m „ = (2ir/x 1 y 2 ){my 2 ; -mx 2 + nxi) 



(9) 
(10) 



(m, n = 0, ±1,±2, ...; triangular lattice: x\ = a, 
x 2 = xi/2, y 2 = xiv / 3/2; square lattice: x\ = y 2 = 
a, x 2 = 0). In (8) Q,A{x,y) is the supervelocity of 
the Abrikosov B c2 solution, which satisfies 



V x Qa = 



B-$o^<5 2 (r-R) 



R 



(11) 



where S 2 (r) = S(x)5(y) is the 2D delta function. 
This shows that Qa is the velocity field of a lat- 
tice of ideal vortex lines but with zero average rota- 
tion. Close to each vortex center one has QaM ~ 
r' x z/(2nr' 2 ) and w(r) oc r' 2 with r' = r R. We 
take Qa from the exact relation 



Qa(t) = 



Vuja x z 
2 koja 



(12) 



where u)A{x,y) is the Abrikosov B c2 solution given 
by a rapidly converging series (6) with coefficients 
a K = -(-) mn + m+n cxp(-K^ n x iy2 /8TT). The solu- 
tions oj(r) and B(r) are then computed by iterat- 
ing equations for the coefficients sk and 6k that 
are derived from the GL equations 5F/5u> = and 
5F/SQ = 0. This method and the obtained GL so- 
lutions are presented in detail in [12,13]. 

The field density P(B) of the ideal triangular vor- 
tex lattice is shown for several b = Bj B c2 values in 
Fig. 1 (k = 1/V2) and Fig. 2 (k = 1.4). For larger 
k the P{B) look similar to Fig. 2. For small b and 
large k the P(B) obtained from the London approx- 
imation are depicted in [11]. The maximum and the 

2 equal minima per unit cell yield two steps in P(B) 
at B — B max (=1) and B = B m - ln (=0), and the 

3 equal saddle points yield a logarithmic infinity at 
B = S sad where P(B) a - In \B - B sad \. 

The variance of the magnetic field, a = ([B(x,y) — 
B] 2 ) = X)k#o^K' i s plotted for the entire ranges 
of b and k in Fig. 3. In the low-field range 
0.13/k 2 « d « 1 one has for the triangular lat- 
tice the London limit a = 0.00371$g/A 4 (upper 
frame in Fig. 3), at very small b -C 0.13/k 2 one has 
a = (6k 2 /87t 2 )<I>o/A 4 (dash-dotted straight lines in 
Fig. 3), and near 6=1 one has the Abrikosov limit 
cr = 7.52 • 10- 4 ($ 2 /A 4 )[k 2 (1 - b)/(n 2 - 0.069)] 2 , 
approximately valid even at 6 ^ 0.3, see [13]. Note 
that the usual London limit for a applies only in 
a narrow range of small (but not too small) 6 and 
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for k > 50. The same is true for the London limit 
of the magnetization curve, where the often used 
"logarithmic law valid at B c \ « J}„ « B c2 for 
M{B a ) = B — B a has a small range of validity [13]. 



8 
7 
6 

S" 5 

ST 

4 

3 
2 
1 




7 
6 

m 5 

Q_ 

4 
3 
2 
1 




x 2 '/x 2 = 1 










0.4 


» 0.95 


m 0.3 


/ II 09 


o_ 






0.2 


Wl A - 8 


0.1 




\0.7 o 



i 0.6 



Vortex lattice sheared along x 



x 2 7x 2 = 1,.95, 
k = 0.707, b=0.5 



.8, .7, .5, .3, 0, 



0.05 



0.1 



0.15 



0.2 



0.25 



Fig. 4. Vortex lattice sheared along the x axis. 
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Fig. 5. Vortex lattice sheared along the diagonal x = y. 

3. Sheared vortex lattices 

The above Fourier method applies to any vortex 
lattice symmetry with vortex positions ~R mn (9), 
also to sheared triangular lattices and to square 
and rectangular basic cells. In Fig. 4 the field den- 
sity P(B) is shown for a lattice sheared away from 
the ideal triangular lattice by decreasing in Eq. (9) 
the length x 2 = 0.5 by a factor c = x' 2 /x 2 = 1 
(triangular), 0.95, 0.9, . . ., (rectangular lattice). 



With y 2 = V3/2 the shear strain in these cases is 
7 = (1 — c) £2/2/2 = (1 — c)/\/3. Note that for this 
shear the saddle point peak splits into 2 peaks, i.e., 
there are now two different types of saddle points in 
B{x, y), but still one maximum and 2 equal minima. 
Figure 5 shows a different orientation of shearing the 
triangular lattice, namely, now the lengths y 2 and x\ 
are changed by a factor c = y 2 /y 2 = x\/x\ = 1.06, 
1.03, 1, 0.97, 0.94, 0.91 such that the unit cell area 
x iV2 = x[y' 2 does not change (no compression). This 
corresponds to a shear strain of size 7 = 2(1 — c) ori- 
ented along the diagonal x = y. One can see that for 
c > 1 the saddle point peak splits into 3 peaks (i.e. 
all 3 saddle points now occur at different B) while 
for c < 1 there occur 2 different saddle point peaks. 

One notes that even very small shear of the vortex 
lattice causes pronounced change in the field prob- 
ability P(B). Small shear costs very little energy 
since the shear modulus c^e of the vortex lattice is 
much smaller than its compressional modulus en or 
its tilt modulus C44. One has approximately [3] 



c n (k) 



B 2 dB a 



1 



+ C 6 6 



C66 



^0 dB (l + k 2 \> 2 ) 
BB c2 (1 - b) 2 (2n 2 - 1)2k 2 



8k 2 ho (2k 2 
B 2 1 

C44(fc) = — 

UO 1 



+ 



1 + 1/Pa) 2 
B{B a - B) 



(13) 
(14) 
(15) 



W 1 + fc 2 A' 2 no 

with A' = A/Vl - b. In c 66 , (3 A = 1.1596 is the 
Abrikosov parameter of the triangular lattice (the 
square lattice is unstable and thus has negative c%&), 
and the factor (2k 2 — 1) means the shear stiffness 
of the vortex lattice is zero in superconductors with 
k = 0.707 (pure Nb). 

An interesting property is the dependence of en 
(13) and C44 (15) on the magnitude k = |k| of the 
wave vector k = (k x ,k y ,k z ) of spatially periodic 
strain, which means the elasticity of the vortex lat- 
tice is non-local. In the limit of uniform stress, k — > 
0, these expressions reproduce the known values of 
the uniform compression and tilt moduli obtained 
by thermodynamics, Cn — c 6 6 = (B 2 / ^ )dBa/dB, 
C44 = BB a //j, . However, when the wavelength of 
the periodic compression or tilt decreases, i.e., k 
increases, these moduli decrease. This means, the 
vortex lattice is softer for short-wavelengths com- 
pression and tilt than it is for long wavelengths. In 
anisotropic superconductors these moduli at finite 
wavelengths are even smaller [14] and the vortex lat- 
tice is softer and can be distorted and melted more 
easily in high-T c superconductors. 
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Fig. 6. Contour lines of order parameter u}(x,y), Eq. (17), 
identical to the contours of induction B(x,y), Eq. (16), for a 
vortex lattice with one vacancy. The displacements s„ of the 
relaxing vortices are shown as short bold lines connecting 
two dots. At all vortex positions ui(x, y) has a minimum and 
is zero, but at the origin x = y = 0, ui is maximum since a 
vortex was removed from there. 
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Fig. 7. P(B) Field density P(B) of a vortex lattice with var- 
ious vacancy concentrations 1/N : a) no vacancy, b) TV = 9, 
c) N = 6, d) N = 4, c) N = 3. GL solution near B c2 . 

4. Vortex lattice containing vacancies 

As an example for structural defects I consider 
a vortex lattice (spacing a) with a super lattice of 
vacancies with spacing Na, N = 2, 3, 4, .... This 
problem was solved in [2] both within London theory 
(at 6 <C 1) and GL theory near 6=1 where 



holds, i.e., the shapes of ui(x, y) and B(x, y) are the 
same. For this vacancy lattice near b = 1 one has 



w(r) = ci 



u A (r) -i-r wa[(i" - R„ 



n 



)/N] 



cj a [(v-R v )/N] 



(17) 



where u>A(x,y) is the Abrikosov B C 2 solution given 
below Eq. (12), c\ is a normalization constant, the 
product is over all vortex positions R„ = R mn 
within the super cell, and the vortex displacements 
s„ are chosen such as to minimize the free energy 
and the Abrikosov parameter f3 = (lu 2 )/(u>) 2 > 1. 
This relaxation of the vortex positions around the 
vacancy yields an uj(x, y) with nearly constant spa- 
tial amplitude, i.e, the maximum w(0,0) at the 
vacancy position has about the same height as all 
the maxima of u> between the vortex positions. 

Figure 6 shows the contours of ui (and thus of B) 
for a vortex lattice with one vacancy (limit N — > oo) 
at the origin and with central symmetric displace- 
ments s„ = -R y [v / 3a 2 /(47ri?2) + o.068a 4 /i?4]. Thc 
field density P{B) of vortex lattices with various va- 
cancy concentrations \/N 2 is shown in Fig. 7. The 
new peaks indicate that new saddle points and min- 
ima of B{x,y) (i.e., maxima of u) appear near the 
vacancy, as seen also in Fig. 6. 
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triangular vortex 
lattice in bulk 




Fig. 8. The field density P(B) of thc triangular vortex lattice 
in a superconducting film of thickness d = 1.2a = 6A in a 
perpendicular field B a = B for re = 0.707, b = B/B c 2 = 0.6, 
from GL theory. Shown are P(B) for bulk (2D, solid line) 
and for film (3D, dashed line). The inset shows the magnetic 
field lines in and near the film; the dashed line marks the 
film surface z = d/2 = 0.6a, a = vortex spacing. 



B(r) = B + [(u) — w(r) ] B c2 /(2k 2 



(16) 
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5. 3D vortex lattice in films 

The Fourier method of Sec. 2 may be general- 
ized to superconducting films of arbitrary thickness 
d (\z\ < d/2) containing a vortex lattice that is pe- 
riodic in the (x, y) plane. For this one has to add to 
F, Eq. (5), the energy F stTay of the magnetic stray 
field outside the film and has to use z dependent 
trial functions, e.g., (6)- (8) with Fourier coefficients 
depending on 3D vectors K = (K x , K y , K z ) [15]. 
When the applied field B a is along z (perpendic- 
ular to the film plane) then the short vortex lines 
are straight and along z. The inset of Fig. 8 shows 
how the magnetic field lines of the vortices become 
less dense when they approach the film surface, they 
"mushroom" and go over smoothly into the stray 
field. 

Though in the depicted example (b = 0.6, k = 
0.707, d = 1.2a «6A« 4£) this widening of the 
field lines is only weak, it still has a strong effect on 
the field density P(B) shown as dashed line, namely, 
the jumps of B m ; n and -B max and the saddle-point 
peak are smeared and shifted towards B, and P(B) 
is increased (has a hump) between B = i? sa d and 
B. These features are expected since for very thin 
films with d < A the field amplitude is strongly re- 
duced inside the film and P(B) narrows to a line po- 
sitioned at B — B. This P(B) should be observable 
by /iSR in specimens composed of many thin layers 
separated by a distance > a/4 where the stray-field 
modulation amplitude cx cxp[— 2tt(\z\ — d/2) /a] has 
almost vanished. 

In such infinitely extended films one has B = 
B a since all field lines have to pass the film. The 
magnetization M of the film, therefore, cannot be 
calculated as a difference of fields, but one has to 
take the derivative of the total free energy, M — 
—d(F + Fstray /d)/dB. A more elegant method cal- 
culates M by Doria's virial theorem directly from 
the GL solution for the film, with no need to take 
an energy derivative [16]. 

6. Random perturbations 

In real superconductors randomly positioned 
weak pinning centers, or random pinning forces, 
may lead to more or less random small displace- 
ments of the vortices from their ideal lattice posi- 
tions. As shown in [11], in bulk superconductors 
this leads approximately to a convolution of the 
ideal-lattice P(B) with a Gaussian. Computer sim- 



ulations of this problem based on London theory 
(i.e., pairwise interacting vortex lines and linear 
superposition of vortex fields) are presented in [11]. 
Interestingly, while disorder of a 2D vortex lattice 
broadens P(B) and its singularities, disorder in 
the 3D point-vortex (or pancake-vortex) lattice oc- 
curring in layered high-T c superconductors [6,17] 
typically will lead to narrowing of P(B) [18] since 
the vortex lines (pancake stacks) become wider. A 
further effect that contributes to the broadening 
of P(B) is the (quantum) diffusion of muons after 
they have stopped, e.g., in ultrapure Nb [19-21]. 

Improved pinning simulations using GL theory 
and considering also thermal fluctuations of vor- 
tices are desirable, as well as microscopic calcula- 
tions going beyond the GL approach. From BCS- 
Gor'kov theory it is shown in [11,19] that for pure Nb 
near B c i the P{B) depends on temperature T, and 
at T <C T c it looks quite different from the GL re- 
sult valid near the critical temperature T c since then 
B(x, y) has sharp conical maxima and minima, and 
two saddle points with three-fold symmetry yielding 
an infinity of the form P(B) oc \B — S sa d| _1 ^ 3 - 
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